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^ ' Quantum effects in a system of coupled atomic and molecular Bose-Einstein condensates in the 

, ^ , . framework of a two-mode model are studied numerically and analytically, using the discrete WKB 

approach. In contrast to the mean-field approximation, the WKB analytical results are in a very 
^ good agreement with numerical results. The quantum fluctuations of the atomic and molecular 

■ populations are calculated, and found to be of the same order of magnitude as their mean values. 
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(3 I. INTRODUCTION 



There is a growing interest in super-chemical P| properties of weakly interacting gases in the Bose-Einstein con- 



I ■ densate (BEC) state. One of the major goals in this field is to produce a molecular Bose-Einstein condensate. To 

' . . . . 

^ ■ date, this has been achieved by creating it from an atomic BEC. Two routes are mainly used: ultra-cold molecules 

O ■ are formed by photo-association of atoms in a BEC 0, 0] or they are produced by applying a time- varying magnetic 
^ field near a Feshbach resonance Q, |H 1^ - 

, , , One of the first theoretical attempts to describe the production of a molecular BEC via Raman photo-association 
of an atomic BEC has been put forward in Ref. Rj- There, the possibility of coherent formation of a molecular BEC 
Ci is treated using a parametric field theory. In Ref. 8], a two-mode quantum Hamiltonian was proposed, that accounts 
for the statistics of both the atoms and the molecules. In this model, large amplitude nonlinear oscillations between 
atomic and molecular condensates are observed, that are damped away at long times. Using the same two-mode 
model, the dynamical evolution of an atomic-molecular BEC has been investigated in Ref. Q and compared to the 
mean- field theory predictions. One of the conclusions drawn in this paper is that, large amplitude atomic- molecular 
oscillations, which were predicted in Ref. 8] on the basis of mean-field results, are damped by the rapid growth of 
fiuctuations near the dynamically unstable molecular mode. Within the same framework, the bi-stability and quantum 
fiuctuations in coherent photo-association of a BEC have been considered in Ref. The number statistics of the 
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two-mode model have also been compared to the number statistics of normal atomic Fermi gases and of Fermi systems 
with pair (BCS) correlations Quite recently a coherent superposition of atoms and molecules condensate was 

observed in a Rb Bose-Einstein condensate by using a stimulated Raman adiabatic passage Q| 

In this paper, we show that, in the two-mode model proposed in Ref. the expectation value of the population 
imbalance between atoms and molecules, as well as the level fluctuations, depend drastically on the initial value of the 
imbalance. We study this pro blem numerically and analytically. Our analytical results, based on the discrete WKB 



pro 



method developed in Ref. [13|, are confirmed by the numerical solution of the exact quantum equations. 



II. THE MODEL 

We consider the simplest model of atomic-molecular condensate, consisting of an atomic mode A and a molecular 
mode B. Let b{b^) and a(a^) be the annihilation (creation) operators for the molecular field B and the atomic field 
A, respectively. The two modes are coupled coherently by an association-dissociation process of a Fermi-resonance 
|l4| type, and the effective Hamiltonian has the form 

H = —a'' a+ (b'' aa + a^ aH) (1) 

where the detuning parameter A — 2fia ~ Mm characterizes the difference between the chemical potentials of the 
molecular (/im) and the atomic (/ia) modes. In the case of coherent photo-association of a Bose-Einstein condensate, 
the Fermi-coupling parameter x/v^ (Rabi frequency) depends on the volume Y of the system [isf . 
Using Hamiltonian (J^), the Heisenberg equations for the operators 

h{t)=a\t)a{t), x{t)=h\t)a{tY +a\tfh(t) and yit) = i {h\t) a(tf - a'^ (tf b{t)) , (2) 

are given by 

^mx ^ = -Ay(t), "^^^m-^Amf^^^m-^)- (3) 



dt " dt ' dt 

The total number of atoms in the system, N = a^a + 2b%, is a conserved quantity. 

In the mean-field approximation, (ri^) « (f^)^, where (...) stands for the average over the initial state of the system. 
Neglecting the term x^l which is of order N times less than the remaining terms in the right-hand side of Eq. Q , 
the equations of motion for {x{t), y{t), n(t)) have stationary point (0, 0, 0), which corresponds to the entire population 
being in the molecular mode. However, for < 2%^ i^, where — N/V is the concentration of atoms in the system, 
this state is dynamically unstable j^. It is one of the goals of this paper to go beyond the mean-field approximation 
and to clarify the behaviour of an atomic-molecular Bose-Einstein condensate in the limit of strong Fermi-coupling. 

III. DYNAMICAL EVOLUTION OF THE NUMBER OF ATOMS 



In what follows, we shall use the basis of states 



(6t)»b 



|n„,nfc) - ^^^=^|0). (4) 

Because the total number of atoms, N, is conserved, \na, nb) = \na, {N — na)/2) = \na). Notice that Ua and N have 
same parity. For the sake of simplicity, we assume N even and put N = 2P, = 2p (p e {0, 1, 2, • • • , P}). We look 
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for a solution of the Schrodinger equation H\(j)) = E\(j)) in the form \(f}) = X]^=o ^p\^P)- The coefficients Cp satisfy 
the set of equations 

£;Cp = + % (/(p)Cp_i+/(p+l)Cp+i) , {0,1,2,. •• ,P}, (5) 

V V 

with 

f{p) = ^p(p-i)(P+l-p) . (6) 

The number of possible energies and eigenstates, jt/)'^"))}, is P + 1. In a system containing initially tlq — 2pa 

atoms, the expectation value of the fraction of atoms, z{t) ~ n{t)/N, evolves according to 



z(t) ^{nole'^^ne '"*\n„) = z + ^ /a,a'K)cos w^.c'* , (7) 

where 



N 

a<a' 



represents the constant (dc-) part of z{t) and where 



iV 



are the respective intensities of the P{P + l)/2 frequencies uja^a' — E^"" ' — P^"-*, a' > a, appearing in the system. 



IV. DISCRETE WKB SOLUTION 

To gain some insight into the dynamics of the system, we use the so-called discrete WKB-approach (see Ref. [l^ l 
to solve the three-term recurrence equation ^ in the limit where P 3> 1. We introduce a new variable x = p/P 
and consider the coefficients C{p) = C{xP) = c{x) as functions of the continuous variable x. Now, up to order 2 in 
e = 1/P, Eq. © reads 

(X-Sx) c{x) - (e (d^F) sinh (e d^) + 2F cosh (ed^)) c{x) = (10) 

where E/ (^/2^x-P), S = A/ (V2^x) and F{x) = p-^/"^ f{xP +1/2) w We now use the WKB method 

to treat equation (|10|) . We look for a solution in the form c(a;) = exp || S'(a;)}. Noting that exp {e9j;} exp||S'(x)} = 
(l + i|S"') exp {iS"} exp S'(a;)} + 0{e'^), where ' = dx, and expanding the "action" S{x) as a power series in e, 
S = So + ieSi, we finally obtain 

AC -I 

: cosS*^ = ~ ^ , and : S*! = -- In (P sinS*^) . (11) 

2r 2 

The "classically allowed region" corresponding to this problem is determined by the inequality 

X-{X) <x <x+{X), (12) 



where the turning points x±{\) are solutions to 

{X-5xf ^ 4p2(a;) 



(13) 
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From , the semi-classical quantization rule for the spectrum can be shown to be given by 

S{X)= J dx arccos ^ 2F(1)^ ) ^ (" + ^) ^ ' '^^^^ 

x_(A) 

where a is a non-negative integer |l3| . This relation yields the number of states with energy lower than A, a(A), 
whose derivative is the density of state given by 



p(A) = 



da[\) 


-hi 


dX 





v[x. A) 



where the "velocity" v{x, A) is given by v{x, A) = y AF'^{x) — (A — 5 x)^. An approximate expression for the coefficients 
inside the allowed region is 

X 

cix, A) — — , cos I - 5*0(0;, X) + 9] where 5*0(2:, A) = / dx' arccos I — — ) . (16) 

yJv{x,X) \e J J \2F{x') J 

X-{\) 

The angle 6 depends on the position of the turning points [fill] but its expression is not needed in what follows. We 
will consider that, outside the allowed region, c{x, A) decays fast enough to be neglected. The normalization constant 
A is determined from the relation 



l=.^|Cp(ii;)|2~P J -^cos^ (^-^So{x,X)+ 9^ dx^P^A^piX). 



(17) 



P=0 ._(A) 



To obtain the last equality, we have replaced the rapidly varying cos'^ (5o/e) by its average value 1/2. This yields 
eventually 



c(x. A) = < 



■cos ( i5o(a;. A) -h6i ) , x^{X) < x < x+{X) 



TrPpiX)v{x,X) Ve / ' - ' (^g^ 

, otherwise. 



V. ATOMIC AND MOLECULAR POPULATIONS 

As already mentioned in the introduction, mean-field theory predicts the existence of large-amplitude coherent 
oscillations between atomic and molecular phases while numerical simulations Isl suggest that these oscillations 
are strongly damped. It is therefore important to clarify the long-time behaviour of the atomic-molecular BEC. To 
do so, we propose to evaluate the dc-component of the atomic population, z. It follows from Eqs. and (|18f) . that 
this quantity can be expressed as 

X-tixo) x+W X+(xo) 

" P A_(a;o) X-{\) X^{xo) 

where 

x+iX) 
x_(A) 
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In deriving Eqs. H19|l and H20() . we have taken into account the fast oscillating nature of Cp{E) as a function of p for 
E fixed (see FigQ), and therefore, cos^ Sq) has been systematically replaced by 1/2. Integrations in Eas. lfT3|) and 
(|20|1 can be carried out (see e.g. [iTl) and as a result, we find the dc-component of the atom population to be 

/ ^ xt{\) + {x+{\) - Xt{X)) (21) 



TT y [ K(m) J ?j(a;o,A) 

where xt{X) is the third root of Ea. H13|l {xt < x- < x+), K(m) (E(rn)) is the complete elliptic integral of the first 
(second) kind Q| with parameter m — (a;^ — a::_)/(x+ — xt), and A± = 6 xq ± 2_F(a:o). For weak detuning, \S\ < 
VI — a;o, expression H21|l is well approximated by 



2 



arcsm 



( VA+/(A+-A_)) arcsin ( v/A_/(A_ - A+)) 



TT 



"^^"^ In (8(1 -(52)/^^) + ln(8(l-<52)/|A_|) ' ' ^^^^ 



Fig. m shows the dc-component of the expectation value of the molecular mode population z„ioi = (1 — z) /2, in 
the resonant case ^ = 0, as a function of the initial value of the atomic mode population zq = z{0). The agreement 
between the numerical solution to Eq. jSjl and the analytical result (|21|) obtained within the discrete WKB approach 
is quite remarkable. We also present here the dc-component of {n(t)) obtained from the set of equations Q in the 
mean-field approximation. In the resonant case 6 = the atomic mode population is governed by the mean-field 
equation 

cPz „ , , 

— + 3z2 - 2z = 23 

ClT'' 

where r = X\P^^ is a rescaled time. Under the initial conditions 

dz 

2(0) = 20, ^ =0 

the solution to Eq. H23|l has the form 

z = zo- ZQ^l- q+) sn^ (^zq \J 2o(1 - <l-) t (24) 

where sn(u|/^) is the elliptic sine function with modulus /i = (1 — g-|-)/(l — (/_) and where 

1 

2~o 



q± = t;— ( 1 - Zo ± \/ 1 + 2^0 - 32o 



It is straightforward to obtain from Eq. (|24f) that the dc-component of the atomic population in the mean-field 
approximation has the form 

2™/ = 2o(^g-+ (l-9-)||^). (25) 

Fig. 121 shows that the numerical and analytical WKB-solutions follow closely the mean-field solution (dashed line) 
near the ends of the interval zq G (0, 1) ( for small and large initial atomic populations). However, in the main part 
of the interval, there is a significant discrepancy between numerical solution and mean-field solution. Note that the 
approximate formula for the dc-component given in 122() is also in a rather good agreement with these results: the 
relative error does not exceed 0.08. Fig. |21 displays the dc-component of the molecular population versus 5 and zq. 
As observed from this graph, for a wide range of initial atomic populations, zq, a significant part of the population is 
on average in the molecular mode. 
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VI. POPULATION FLUCTUATIONS 

To clarify the relevance of the molecular population behaviour shown in Figs. [21 and 13 to other physical properties 
of the system, we need to estimate the fluctuations of these quantities. To this end we consider the dc-component of 
the mean square fluctuation of the number of atoms 

1 



{Any ^ {h^{t)) (hit)) , (h^t)) ^ hm - (n^(t)) dt. (26) 

1 — *oo 1 







From Eqs. ^ we find that the equation for the operator of the number of atoms n{t) is 

dMt) • s^n{tf+(A^-2^N)h{t)-^N-2AH = (27) 



dt^ V " ' V V J ^ ' V 

where H is the Hamiltonian of the system Now, averaging this last equation with respect to any quantum state 
and taking its dc-part, yields the following result exact for any quantum state and any total number of particles N 

From Eqs. and (|28|l . we find that, for an initial state given by Eq. I^J with Ua — uq {nt — {N — no)/2), the mean 
square fluctuation of the number of molecules in the atomic-molecular BEC is 

with the dc-component of the atom population z given by Eq. H21|) . The relative fluctuation of the number of 



molecules is given by the expression v = y {Ant)'^ ^ (nb). From Fig. ^ it can be seen that the fluctuations of the 
number of molecules are of the same order of magnitude as their mean values. It is worth recalling that the mean 
square fluctuations of the number of particles in an ideal gas is of the order of and are exceedingly small in 
macroscopic assemblies. On the other hand, the mean square fluctuations of the populations of individual energy 
states in an ideal Bose-Einstein gas are known to be greater than 1 (see e.g. ^3). The two-mode atom-molecule BEC 
is then an intermediate system with mean square fluctuations both less than 1 (see Fig. ^ and greater than 
Given the large fluctuations experienced by the populations around their average values, there is, strictly speaking, 
no convergence (in time) towards any state. 



VII. CONCLUSION 



In summary, we have considered the two-mode model of atom-molecule Bose-Einstein condensate, and studied 
quantum effects on atom-molecule population oscillations. We used numerical simulations and analytical approaches. 
We found numerically eigenvalues and eigenvectors of the full N-body problem. Analytical calculations were performed 
by using the discrete WKB approach. We found excellent agreement between the results of direct numerical simulations 
and analytical results. In contrast, the results of the widely used mean-field approach deviate significantly from the 
numerics for a wide range of atom concentrations. We studied also quantum fluctuations in coherent atom-molecule 
transformations. We found that that the fluctuations of the number of atoms (molecules) are of the same order of 
magnitude as their mean values. Thus, the two-mode BEC is seen to be in an intermediate regime: molecules and 
atoms can still be treated as separate phases, but their populations are highly fluctuating. 
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FIG. 1: Coefficients Cp{E) (arbitrary units) obtained from Eq. ll^J for 5 = and = 500 (P — 250) and placed at the 
corresponding energy E ~ 1580.8 (thin dashed line). The allowed region is enclosed by thick upper and lower lines. Coefficients 
Cp{E) are clearly oscillating in this region and rapidly decaying outside. 
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FIG. 2: Normalized dc-component of the molecular population Zmoi ~ {l~z)/2 versus the normalized initial number of atoms, 
z(0) for a zero detuning parameter, S — 0. Exact numerical results are obtained from Eqs. and ^ for A'' = 500 cannot 
be distinguished from the analytical expression 12111 shown as the solid line. The dashed curve represents the same quantity 
obtained in the framework of the mean field approach 
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FIG. 3: Dc-component of the molecular population, Zmoi = (1 — a) /2, versus the detuning parameter S and the initial value 
of the atomic population zo- 
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FIG. 4: Dc-component of the relative fluctuation of the number of molecules, v — \J (Arat,)^^ (uj), for different values of the 
detuning parameter 5: 5 = Q (solid line), 5 = 0.5 (dashed line) and 6 = 0.75 (dot-dashed line) . 



